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ABSTRACT 


Stiffened  plates  are  designed  for  minimum  volume  subject 
to  constraints  on:  Vor  Mises  maximum  stresses,  nodal  dis¬ 
placements,  height  to  thickness  ratios  of  frame  webs,  and 
width  to  thickness  ratios  of  frame  flanges-  Design  variables 
are  plate  thicknesses  and  stiffener  dimensions. 

A  finite  element  analysis  program  is  developed  for  the 
design  of  stiffened  plates  using  numerical  optimization 
techniques.  The  program  may  be  used  as  a  stand  alone  analy¬ 
sis  tool  or  may  be  coupled  to  an  optimizer  of  user's  choice. 
Rectangular  plate  elements  and  frame  elements  are  used  for 
the  idealization  of  stiffened  plates. 

Design  examples  are  presented  to  demonstrate  the  design 
method. 


4 


^  *■*  -  *r‘  ^ 


TABLE  OF  CONTENTS 


I.  INTRODUCTION -  10 

II.  FINITE  ELEMENT  METHOD  -  13 

A.  GENERAL -  13 

B.  FINITE  ELEMENT  DISPLACEMENT  APPROACH  -  14 

1.  Element  Analysis  -  14 

2.  Direct  Stiffness  Method  -  17 

C.  FINITE  ELEMENTS  USED -  17 

1.  Rectangular  Plate  Element  -  18 

2.  Frame  Element -  18 

D.  OFFSETTING  OF  RIGID  LINKS  -  24 

III.  OPTIMIZATION  TECHNIQUES  -  28 

A.  INTRODUCTION -  2  8 

B.  DEFINITIONS -  28 

C.  COPES/CONMIN -  31 

D.  DESIGN  PARAMETERS -  32 

1.  Design  Variables  -  32 

2.  Objective  Function  -  32 

3.  Constraints -  33 

IV.  DESIGN  EXAMPLES -  35 

A.  DESIGN  CASE  I— THE  UNSTIFFENED  PLATE -  35 

1.  Case  IA  with  A  Concentrated  Load -  35 

2.  Case  IB  with  Uniform  Distributed  Load  -  38 

B.  DESIGN  CASE  II— THE  STIFFENED  PLATE -  38 


5 


40 


1.  Case  I IA--Rect angular  Type  Frame 

Stiffener  - 

2.  Case  IIB — T  Type  Frame  Stiffener -  43 

C.  DESIGN  CASE  III— THE  HATCH  COVER -  4  7 

1.  Case  IIIA — Rectangular  Type  Frame 

Stiffener -  47 

2.  Case  I  IIB — T  Type  Frame  Stiffener -  52 

V.  CONCLUSIONS  AND  RECOMMENDATIONS  -  57 

A.  CONCLUSIONS -  57 

B.  RECOMMENDATIONS -  58 

LIST  OF  REFERENCES  -  59 

INITIAL  DISTRIBUTION  LIST  -  60 


6 


*'  ■ 


LIST  OF  TABLES 


I.  Summary  of  Design  Case  I -  3  8 

II.  Design  Case  IIA — Design  Variables  and 

Nodal  Connectivity -  41 

III.  Summary  of  Design  Case  IIA -  42 

IV.  Design  Case  IIB — Design  Variables  and 

Nodal  Connectivity -  44 

V.  Summary  of  Design  Case  IIB -  45 

VI.  Design  Case  IIIA — Design  Variables  and 

Nodal  Connectivity -  49 

VII.  Summary  of  Design  Case  IIIA -  31 

VIII.  Design  Case  IIIB — Design  Variables  and 

Nodal  Connectivity -  53 

IX.  Summary  of  Design  Case  IIIB -  56 


7 


T 


LIST  OF  FIGURES 


1.1  Stiffened  Plate  Model  -  11 

2.1  Rectangular  Plate  Element  -  19 

2.2  Frame  Element -  21 

2.3  Stiffener  Options  -  22 

2.4  Rigid  Links -  26 

4.1  Design  Case  I--The  Unstiffened  Plate  Model  -  36 

4.2  Design  Case  I-~Design  Variables  -  37 

4.3  Design  Case  II — The  Stiffened  Plate  Model  -  39 

4.4  Design  Case  III--The  Hatch  Cover  Model  -  48 


3 


ACKNOWLEDGMENT 


I  would  like  to  express  my  sincere  appreciation  to 
Professor  G.  N.  Vanderplaats  for  his  endless  patience  and 
willingness  to  stop  whatever  he  was  doing  and  provide  much 
needed  guidance.  His  positive  attitude  and  encouragement 
were  a  tremendous  support. 

I  would  also  like  to  thank  Professor  Y.  S.  Shin  for  his 
most  efficient  assistance  and  the  faculty  of  the  Mechanical 
Engineering  Department  for  the  instruction  I  received  during 
my  course  of  study  at  the  Naval  Postgraduate  School. 

Finally,  I  would  like  to  thank  my  loving  wife,  Hyun  Sun, 
and  my  cute  daughter,  Seo  Young,  for  their  faith,  encourage¬ 
ment,  and  prayers  during  the  past  year.  They  provided  much 
needed  motivation  to  keep  on  plugging  when  things  got  sticky. 


INTRODUCTION 


I . 

The  development  of  high-speed  digital  computers  has  made 
possible  significant  changes  in  the  structure  design  process. 
One  of  these  is  the  availability  of  various  mathematical 
programming  metnods  for  use  in  design  optimization.  The 
computer's  speed  allows  the  designer  to  now  consider  a  much 
wider  range  of  design  alternatives.  The  optimization  proce¬ 
dure  provides  a  means  of  systematically  choosing  from  among 
these  alternatives  based  on  some  predetermined  rational 
criterion. 

Even  when  the  selected  numerical  method  is  able  to  ar:  re 
at  the  optimum  design,  the  result  is  only  as  good  as  the 
design  model.  Here  there  is  an  even  greater  need  for  design 
experience  and  sound  engineering  judgment.  The  design  model 
must  be  carefully  developed  to  realistically  represent  the 
design  in  question. 

The  finite  element  model  for  stiffened  plate  is  shown  in 
Fig.  1.1,  where  the  plate  has  been  descritized  by  rectangular 
plate  elements  and  the  stiffener  by  frame  type  elements.  The 
eccentricity  of  stiffener  is  transformed  to  the  linked  nodal 
point  by  applying  a  linear  equation  that  realtes  displacement 
degrees  of  freedom. 

The  purpose  of  this  thesis  is  to  develop  a  finite  element 
analysis  program  for  stiffened  plates,  and  to  design  the 
optimum  stiffened  plate  by  coupling  two  programs;  the 
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analysis  program  and  ageneral  purpose  non-linear  optimizer, 
COPES/CONMIN  [Refs.  1-2]. 

The  remaining  chapters  of  this  thesis  are  outlined  as 
follows : 

Chapter  II  presents  briefly  the  finite  element  method 
used  in  the  analysis  program  for  stiffened  plates. 

Chapter  III  presents  the  basic  concepts  cf  the  optimi¬ 
zation  methods  used  in  the  COPES/CONMIN. 

Chapter  IV  presents  design  examples. 

Chapter  V  offers  conclusions  and  recommendations. 
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II.  FINITE  ELEMENT  METHOD 


A.  GENERAL 

The  finite  element  method  is  now  well  established  as  an 
engineering  tool  of  wide  application.  The  fundamental 
concept  of  the  finite  element  method  is  that  any  continuous 
quantity,  such  as  pressure  or  displacements,  can  be  approxi¬ 
mated  by  a  discrete  model  composed  of  a  set  of  piecewise 
continuous  functions  defined  over  a  finite  number  of  sub- 
domains.  The  piecewise  continuous  functions  are  defined 
using  the  values  of  the  continuous  quantity  at  a  finite 
number  of  pcrnrs  in  its  domain. 

The  formulation  of  the  finite  element  method  can  be  traced 
to  energy  procedures,  principally  the  minimum  potential  energy 
method  and  the  minimum  complimentary  energy  method.  The 
minimum  potential  energy  method  is  associated  with  assumed 
displacement  parameters  as  unknowns  and  is  usually  termed 
the  "displacement"  or  "stiffness"  method.  On  the  other  hand, 
the  minimum  complimentary  energy  method  dealt  with  stress 
parameters  and  is  termed  the  "flexibility"  or  "force"  approach. 
The  ease  with  *rtiich  a  continuous  displacement  pattern  can  be 
prescribed  (compared  to  the  alternative  approach  of  forming 
an  equilibrating  internal  force  field)  has  aided  the  wide¬ 
spread  use  and  development  of  the  finite  element  displace¬ 
ment  approach.  The  displacement  model  and  the  stiffness 
analysis  are  employed  in  the  analysis  program  developed  here. 
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This  chapter  will  briefly  present  some  of  the  general 
concepts  of  the  finite  element  method  used  in  the  analysis 
program. 

B.  FINITE  ELEMENT  DISPLACEMENT  APPROACH 

The  displacement  formulation  involves  derivation  of  the 
stiffness  matrix  of  each  individual  element.  The  stiffness 
matrix  of  the  entire  assembled  structure  is  then  obtained  by 
super-position.  This  matrix,  along  with  the  prescribed 
displacement  boundary  conditions  and  loads,  is  used  for  the 
solution  of  displacements  and  stresses. 

1 .  Element  Analysis 

For  the  structural  applications  at  least,  the 
governing  equilibrium  equations  can  be  obtained  by  mini¬ 
mizing  the  total  potential  energy  of  the  system.  The  total 
potential  energy,  t  can  be  expressed  as 

TT  =  if  JT  e  dv  -  /  UT  f  dv  -  /  UT  q  dS  (2.1) 

V  '  ~  V''  s' 

where  a  and  s  are  the  stress  and  strain  vectors  respectively, 
u  the  displacements  at  any  point,  f  the  body  forces  per  unit 
volume  and  q  the  applied  surface  tractions.  Integrations  are 
taken  over  the  volume  v  of  the  structure  and  loaded  surface 
area  3. 

The  first  term  on  the  right  hand  side  of  Eq.  (2.1) 
represents  the  internal  strain  energy  and  the  second  and 


14 


third  terms  are  respectively  the  work  contributions  of  the 
body  forces  and  distributed  surface  loads. 

In  the  finite  element  displacement  method,  the  basic 
steps  for  derivation  of  the  element  stiffness  matrix  are: 

a.  Express  the  displacements  to  have  unknown  values  only 
at  the  nodal  points,  so  that  the  variation  within  any 
element  is  described  in  terms  of  the  nodal  values  by 
means  of  interpolation  functions.  Thus 

u  =  N  ue  (2.2) 

where  N  is  the  set  of  interpolation  functions  termed 
the  shape  functions  and  ue  is  the  vector  of  nodal 
displacements  of  the  element. 

b.  Express  the  strains  within  the  element  from  the 
element  nodal  displacements  as 

e  =  B  ue  (2.3) 

where  E  is  the  strain-displacement  matrix  generally 
composed  of  derivatives  of  the  shape  functions. 

c.  Express  the  stresses  relating  to  the  strains  by  use  of 
an  elasticity  matrix  D  containing  the  appropriate 
material  properties,  as  follows: 

a  =  D  s  (2.4) 

d.  Establish  the  equilibrium  equation  of  element. 


Provided  that  the  element  shape  functions  have  been  chosen 
so  that  no  singularities  exist  in  the  integrands  of  the 
functional,  the  total  potential  energy  of  the  continuum  will 
be  the  sum  of  the  energy  contributions  of  the  individual 
elements.  Thus 


=  y 


(2.5) 


where  ue  represents  the  total  potential  energy  of  element  e 
which,  on  use  of  Eq.  (2.1),  can  be  written 


->  / 


ueT  B  u  B  u6  dV  -  /  u6T  NA  f  dV 

...  o  >  .  e- 


V  . 


V 

e 

,eT  .,T 


-  /  uel  n  q  dS 


(2.6) 


where  V  is  the  element  volume  and  S  the  loaded  element 
e  e 

surface  area.  Performance  of  the  minimization  for  the 
element  e  with  respect  to  the  nodal  displacements  uefor  the 
element  results  in 


Sit 

Su 


=  /  (BT  DB)  ue  dV  -  /  NT  f  d  V  -/  NTqdS 


.  e  e  e 
k  u  -  p 


(2.7) 


where 
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f 


Pe  =  /  NT  f  dV  +  /  NT  q  ds 


are  the  equivalent  nodal  forces  for  the  element,  and 


(2.8) 


ke  =  /  Bf  D  B  dV 


(2.9) 


V. 


is  termed  the  element  stiffness  matrix. 

2 .  Direct  Stiffness  Method 

The  real  elastic  structure  is  now  represented  by  a 
finite  number  of  small,  discrete  elements.  Once  their 
approximate  behaviors,  identified  by  their  individual  stiff¬ 
ness  matrices  keof  Eq.  (2.9),  have  been  established,  the 
stiffness  matrix  K  for  the  complete  structure  is  obtained  by 
the  proper  summation  of  each  element  stiffness  matrix  in  the 
structure.  The  summation  of  the  terms  in  Eq.  (2.7)  over  all 
the  elements,  when  equated  to  zero,  results  in  a  system  of 
equilibrium  equations  for  the  complete  continuum.  This 
assembly  process  is  known  as  the  Direct  Stiffness  Method. 
These  equations  are  then  solved  by  any  standard  technique  to 
yield  the  nodal  displacements.  Note  that  K  is  symmetric  and 
positive-definite . 

C.  FINITE  ELEMENTS  USED 

The  finite  elements  used  in  the  analysis  program  will  be 
described  briefly  in  this  section. 
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1 .  Rectangular  Plate  Element 

The  rectangular  plate  element  used  here  is  illus¬ 
trated  in  Fig.  2.1,  where  each  nodal  point  has  6  degrees  of 
freedom,  3  transitional  displacements  and  3  rotational 
displacements.  This  element  with  4  corner  nodal  points  has 
the  element  stiffness  matrix  of  order  24  x  24.  The  corres¬ 
ponding  displacements  of  each  node  u^ , , . . . , u24  will  be 
taken  to  be  positive  in  positive  directions  of  the  xyz- 
coordinates . 

The  main  assumptions  in  the  method  are  that  displace¬ 
ments  are  small  compared  with  plate  thickness,  the  stress 
normal  to  the  midsurface  of  the  plate  is  negligible,  and 
normals  to  the  midsurface  before  deformation  remain  straight 
but  not  necessarily  normal  to  the  midsurface  after  deformation. 

The  assumed  displacement  functions  will  be  taken  to 

be  linearly  varying  in  the  plane  of  the  element.  These 

displacement  functions  will  ensure  both  deflection  and  slope 

compatibility  of  the  adjacent  elements.  The  stiffness  matrix 

of  the  rectangular  plate  element,  which  is  matrix  product 
T 

B  DB  of  Eq.  (2.9)  integrated  over  its  volume,  is  summarized 
by  Przemieniecki  [Ref.  3]. 

Note  that  element  stiffness  matrices  are  formed 
directly  in  the  global  coordinate  system  so  that  no  trans¬ 
formations  from  local  to  global  coordinates  are  required. 

2 .  Frame  Element 

The  frame  element  as  a  stiffener  has  6  degrees  of 


freedom  for  each  node,  such  as  those  of  plate  elements.  The 
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basic  assumptions  are  that  the  stiffener  element  is  straight 
with  uniform  cross  section  capable  of  resisting  axial  forces, 
bending  moments  about  the  two  principal  axes  in  the  plane  of 
its  cross-section,  and  twisting  moments  about  centroidal  axis, 
and  that  its  deflection  due  to  shearing  strains  are  neglected. 
The  stiffener  element  has  2  kinds  of  options  which  are  illus¬ 
trated  in  Fig.  2.3.  The  width  and  thickness  of  flange  of 
the  rectangular  stiffener  element  may  be  referred  as  zeros. 

In  order  to  determine  the  stiffness  property  of  a 
complete  structure,  a  global  coordinate  must  be  established 
for  all  unassembled  structural  elements  so  that  all  the 
displacements  and  their  corresponding  forces  will  be  referred 
to  this  system.  Since  the  element  stiffness  matrices  k  are 
initially  calculated  in  local  coordinates,  suitably  oriented 
to  minimize  the  computing  effort,  it  is  necessary  to  intro¬ 
duce  transformation  matrices  changing  the  frame  of  reference 
from  a  local  to  a  global  coordinate  system.  The  first  step 
in  deriving  such  a  transformation  is  to  obtain  a  matrix  rela¬ 
tionship  between  the  element  displacements  u  in  the  local 
system  and  the  element  displacements  u'  in  the  global  system. 
This  relationship  is  expressed  by  the  matrix  equation 

u  =  T  u'  (2.10) 

where  T  is  a  matrix  of  coefficients  obtained  by  resolving 
global  displacements  in  the  directions  of  local  coordinates. 
The  transformation  matrix  T  is  given  by 
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Figure  2.2  Frame  Element 


Rectangular  type  frame  stiffener 


T  type  frame  stiffener 


Figure  2.3  Stiffener  Options 
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represents  matrices  of  direction  cosines  for  local  x,  y  and 
z  directions,  respectively,  measured  in  global  system  x,  y 
and  z.  Similarly  the  relationship  for  forces  is  expressed 
by  the  matrix  equation 


P  =  T  p' 


(2.11) 


where  p  is  the  force  vector  in  the  local  system  and  p'  is 

the  force  vector  in  the  global  system.  Matrix  A  is  orthogonal; 
-1  T 

that  is,  T  =  T  .  Therefore 


u'  =  T  u 


(2.12) 


?'  =  T  p 


(2.13) 


i 


Let  u'  and  u  be  two  ways  to  describe  the  same  virtual  dis¬ 
placement.  Virtual  work  is 


>Su  ,Tp  ’ 


=  i 


o  u  p 


T  T 

i  u '  T  p 


(2.14) 


So 


5u,T(p'-TTp)  =  0,  p*  =  TTp  (2.15) 

To  transform  the  stiffness  matrix,  we  start  with  k'u'  =  p' 
and  substitute  from  the  preceding  equations. 

k 1 u 1  =  p'  =  TTp  =  TTk  u  =  TTk  Tu'  (2.16) 

Since  this  relation  is  presumed  valid  for  any  u',  we  conclude 
that  the  required  stiffness  transformation  is 

k’  =  TT  k  T  (2.17) 


D.  OFFSETTING  OF  RIGID  LINKS 

One  of  the  most  important  advantages  of  the  finite  ele¬ 
ment  technique  is  that  an  assembly  of  different  structure 
elements  such  as  plates  and  frames  can  be  dealt  within  a 
single  coordinate  analysis.  Usually  the  neutral  surfaces  of 
plate  and  stiffener  are  not  coincident:  the  stiffener  is  on 
one  side  of  the  plate.  A  standard  preliminary  treatment  is 
to  connect  adjacent  plate  and  stiffener  nodes  by  a  rigid 
link,  so  that  degrees  of  freedom  of  the  stiffener  are  replaced 
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by  degrees  of  freedom  of  the  plate.  The  usual  assembly  is 
then  possible.  The  necessary  transformation  is  now  described. 

The  stiffener  element  in  Fig.  2.4  has  the  usual  12  degrees 
of  freedom--6  at  node  A  and  6  at  node  B.  With  reference  to 
these  degrees  of  freedom,  element  load  and  stiffness  matrices 
are  p  and  k.  Similar  degrees  of  freedom  are  used  at  nodes  1 
and  2  of  plate  element  of  rigid  links  A-l  and  B-2.  The 
"master"  degrees  of  freedom  at  node  1  and  "slave”  degrees 
of  freedom  at  node  A  have  the  relation 
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A  similar  expression  is  written  for  link  B-2  by  replacing 
subscripts  A  and  1  by  B  and  2.  The  transformation  vectors 
p  and  matrix  k,  associated  with  degrees  of  freedom  at  nodes  1 
and  2,  are 

(2.19 

(2.20 


The  foregoing  transformation  makes  the  translational 
displacements  depend  on  the  rotational  displacements  ar.d  sc 
introduces  an  unwanted  quadratic  field  into  the  translation al 
displacements.  It  is  found  that  a  frame-stiffened  plate  is 
overly  flexible  in  a  coarse  mesh.  Mesh  refinement  nelps: 
error  is  reduced  by  a  factor  of  four  if  the  number  of  elements 
is  doubled.  The  error  can  be  eliminated  through  the  addition 
of  one  more  degree  of  freedom.  Further  details  of  the  error 
reduction  are  described  in  [Ref.  4], 


P  =  T  p 


T  — 

k  =  T  k  T 


where 


:12*12 


“1 


0  't 


--2 
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III.  OPTIMIZATION  TECHNIQUES 


A.  INTRODUCTION 

In  this  chapter  some  of  the  fundamental  ideas  and  formu¬ 
lation  methods  of  mathematical  programming  are  introduced 
to  understand  the  optimization  method  and  the  optimizer 
COPES/CONMIN.  Design  parameters  used  here  are  then  dis¬ 
cussed.  Fox  [Ref.  5]  and  Hinmelblau  [Ref.  6]  provide  an 
extensive  discussion  of  numerical  optimization  techniques 
and  their  application  to  engineering  design. 

B.  DEFINITIONS 

In  discussing  the  optimization  methods,  the  following 
definitions  will  be  useful: 

Design  variables — the  design  variables  are  the  numerical 
parameters  for  which  values  are  to  be  chosen  in  producing  a 
design.  In  a  structural  problem,  they  might  be  plate  thick¬ 
ness,  frame  dimensions,  ate. 

Objective  f unction--the  objective  function,  is  the  single 
valued  function  with  respect  to  which  the  design  is  optimized. 
In  a  structural  design  problem,  it  might  be  the  weight,  volume 
or  fabricated  cost  of  the  structure.  The  selection  of  an 
objective  function  can  be  one  of  the  most  important  decisions 
in  the  optimum  design  process. 

Constraints — practical  design  problems  are  usually  sub¬ 
ject  to  a  series  of  constraints  which  must  be  satisfied  in 
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order  to  produce  an  acceptable  design.  These  constraints 
may  be  linear  or  non-linear.  If  a  parameter  is  beyond  the 
value  of  a  specified  value,  it  is  said  to  be  violated. 

Side  constraint — the  side  constraint  is  a  constraint 
which  restricts  the  specified  range  of  a  design  variable 
for  reasons  other  than  the  direct  consideration  of  performance. 

Feasible  design — the  feasible  design  (or  acceptable 
design)  is  a  design  which  satisfies  the  specified  constraints. 

Infeasible  design--the  infeasible  design  (or  unaacep- 
table  design)  is  a  design  which  does  not  satisfy  the  speci¬ 
fied  constraints. 

The  general  nonlinear  optimization  problem  is  expressed 
mathematically  as  follows: 

Minimize  F(x) 

Subject  to: 

gj  (x)  <_  0  for  j  =  l,m 

o  n 

x.  <  x.  <  x,  for  i  =  l,n 

1  —  x  —  i 

where  the  vector,  x,  is  the  vector  of  n  design  variables. 

The  objective  function,  F(x),  given  by  Eq.  (3.1),  as  well 
as  the  constraint  functions  given  by  Eq.  (3.2),  may  be  linear 
or  non-linear  functions  of  the  design  variables.  They  may  be 
explicit  or  implicit  functions  of  x,  but  must  have  continuous 
first  derivatives.  If  it  is  desired  to  maximize  F(x) ,  the 
minimization  of  the  objective  function  is  used  since  maximum 


(3.1) 

(3.2) 

(3.3) 
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F(x)  can  be  treated  minimum  of  -F(x) .  The  function  g.. (x) 

~  J 

is  the  set  of  inequality  constraints  to  be  met.  The  m  in¬ 
equality  constraints  must  be  satisfied  to  be  a  feasible 

ll 

design.  Side  constraints,  x^  and  x^,  are  the  lower  and  upper 
limitations  placed  on  the  design  variables.  Side  constraints 
could  be  included  in  Eq.  (3.2),  but  are  treated  separately 
for  efficiency. 

The  n-dimensional  space  spanned  by  the  design  variables 
x^  is  referred  to  as  the  design  space.  As  stated  previously, 
any  design  which  satisfies  the  inequalities  of  Eq.  (3.2)  is 
referred  to  as  a  feasible  design.  If  the  design  violates  one 
more  of  the  inequalities,  it  is  said  to  be  infeasible.  The 
minimum  feasible  design  is  said  to  be  optimal. 

Most  nonlinear  optimization  programs  update  the  vector 
of  design  variables  by  the  iterative  relationship 


xq+1  =  Xq  +  a*  S 


(3.4) 


where  q  is  the  iterative  number,  vector  S  is  the  direction 
of  search  in  the  design  space,  and  a*  is  a  scalar  which 
defines  the  distance  of  travel  in  the  direction  S  during  the 
qth  iteration.  An  initial  design  defined  by  X°  must  be 
defined  and  may  be  a  feasible  or  infeasible  design. 

The  optimization  process  then  proceeds  in  two  steps.  The 
first  is  the  finding  of  S  which  will  improve  the  design 
without  violating  constraints  and  the  second  is  the  determina¬ 
tion  of  a*  which  will  improve  the  design  as  much  as  possible 
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when  moving  in  this  direction.  The  process  is  then  repeated 
until  there  is  no  further  design  improvement,  indicating 
that  this  is  the  optimum  attainable  design. 

C.  COPES/CONMIN 

The  COPES/CONMIN  optimization  program  is  a  general  pur¬ 
pose,  non-linear  optimizer  capable  of  handling  large,  con¬ 
strained  problems.  This  includes  the  conjugate  direction 
method  of  Fletcher  and  Reeves  [Ref.  7]  for  unconstrained 
function  minimization  and  a  modification  of  Zoutendijk's 
Method  [Ref.  8]  of  Feasible  Directions  for  constrained  func¬ 
tion  minimization.  It  has  been  successfully  used  in  connec¬ 
tion  with  structural  optimization  [Ref.  9],  airfoil  design  \ 
[Ref.  10],  aircraft  synthesis  [Ref.  11],  and  numerous  other 
engineering  applications. 

It  was  necessary  to  develop  a  subroutine,  ANALIZ,  which 
for  a  given  design,  would  analyze  a  stiffened  plate,  and  which 
would  be  suitable  for  coupling  with  the  optimizer.  The 
common  block  GLOBCM  is  required  to  couple  the  analysis  sub¬ 
routine  directly  to  the  COPES/CONMIN.  All  variables,  which 
are  used  as  objective  function,  constraints  and  design  varia¬ 
bles,  must  be  listed  in  the  common  statement  and  the  statement 
must  appear  in  each  subroutine  the  variables  are  used  in. 

It  is  used  by  the  optimizer  as  a  catalog  to  identify  where 
the  design  variables,  objective  function  and  constraints  are, 
and  what  purpose  they  fulfill. 
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In  order  to  execute  the  COPES/CONMIN  program  it  is 


necessary  to  provide  formatted  data  for  C OPES/CONMIN,  followed 
by  data  for  the  ANALIZ  program.  Further  details  for  the 
COPES  data  are  explained  in  [Ref.  1]. 

D.  DESIGN  PARAMETERS 

The  design  variables,  constraints  and  objective  function 
used  in  the  design  process  are  discussed  in  this  section. 

1 .  Design  Variables 

The  design  variables  are  plate  thicknesses,  t.,  and 
stiffener  dimensions,  being  the  web  heights  h^,  the  web 
thicknesses  t^,  the  flange  widths  W^,  and  the  flange  thick¬ 
nesses  t-^,  i  =  l,ns,  where  ns  is  the  number  of  the  different 
stiffener  dimension  sets.  Thus,  the  total  number  of  design 
variables  depends  on  the  sets  of  plate  thicknesses  and  the 
sets  of  different  stiffener  dimensions. 

2 .  Objective  Function 

Total  structure  volume  is  considered  as  the  objec¬ 
tive  function  in  the  design  process. 

m  ns 

volume  =  y  v  .  +  y  N .  V  . 
i=l  P1  j=l  ^ 

where  V  ^  is  the  volume  of  the  ith  plate  element,  m  is  the 
number  of  plate  elements,  is  the  volume  of  the  jth  stiffener 

set,  and  Nj is  the  number  of  stiffeners  in  this.  The  number 
of  stiffener  sets  is  ns. 


3.  Constraints 


Design  constraints  are  Von  Mises  maximum  stress, 
nodal  displacements,  height  to  thickness  ratios  of  frame 
webs,  and  width  to  thickness  ratios  of  frame  flanges. 
Stress: 


u  /  i 

a  .  ./a  -  1 

xj  max 


i  =  l,ne 
j  =  l^nj 


Z  ,  , 

a i j^0  max 


i  =  l,ne 
j  =  1/nj 


where  and  are  respectively  the  upper  and  lower  sur¬ 

face  stress  at  the  node  j  of  the  element  i,  and  o  is  the 

J  max 

maximum  allowable  stress.  The  number  of  elements  is  ne  and 
nj  is  the  number  of  nodal  points. 

Displacements : 


u.  ./u 
x]  max 


1  <  0 


i  =  1 ,  n  j 
j  =  l,nc 


where  u^  is  the  displacement  at  node  i  in  the  coordinate 
direction  j  and  umax  is  the  maximum  allowable  displacement 
in  coordinate  direction  j .  The  number  of  coordinates  is 
nc . 

Height  to  thickness  ratio  of  frame  webs: 


h./t  ■  -  5  <  0 

r  Wi  — 


1  -  h./t...  <  0 

l  Wi  - 


i  =  l,ns 
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where  and  are  respectively  the  web  height  and  web 

thickness  of  ith  stiffener  set. 

Width  to  thickness  ratios;  of  flanges: 


-  5  s.  0 


i  =  l,ns 


1  -  w . /t ,  .  <  0 

v  f  i  - 


where  w^  and  t^  are  respectively  the  width  and  the  thickness 
of  flange  of  the  ith  stiffener  set. 


"Vs  *  &*.**+■ 


IV.  DESIGN  EXAMPLES 


This  chapter  will  present  design  examples  of  unstiffened 

and  stiffened  plates.  All  calculations  were  carried  out  on 

the  IBM  3033  model  370  computer  at  the  Naval  Postgraduate 

School.  Material  constants  used  for  all  examples  are  Young's 

£ 

modulus,  E  =  3*10  psi,  and  Poisson's  ratio,  v  =  0.3.  Due 
to  symmetry  only  a  quarter  of  the  plate  was  modelled  in  each 
example . 

A.  DESIGN  CASE  I — THE  UNSTIFFENED  PLATE 

A  quarter  of  the  unstiffened  plate  simply  supported,  as 
shown  in  Fig.  4.1,  was  modelled  using  a  4  *4  mesh  of  the 
plate  elements.  There  are  10  design  variables  as  illustrated 
in  Fig.  4.2.  The  volume  of  the  plate  was  minimized  subject 
to  a  total  of  378  constraints  being  stress  and  nodal  deflec¬ 
tion  constraints,  with  a  =  20,000  psi  and  u  =  1/2  t. 

max  ^  max 

1 .  Case  IA  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  was  applied  at  the 
center  of  the  whole  plate.  The  final  optimization  results 
show  that  there  are  3  critical  deflection  constraints  of 
the  nodes  3,  4  and  10,  and  that  the  plate  thicknesses  of  the 
diagonal  elements  are  much  thicker  than  those  of  off-diagonal 
elements.  The  final  volume  is  363.916  in'*.  The  results  of 
this  case  are  summarized  in  Table  I  for  comparison  with  the 
other  cases. 
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t  (8) 


t  (5) 


tCl) 


TABLE  I 


Summary  of  Design  Case 

I 

Variable 

Concentrated 

Distributed 

Load 

Load 

t(l) 

1.  302 

0.352 

t  (2) 

1.031 

0.263 

t  (3) 

0.573 

0.589 

t  (4) 

1.059 

0.799 

t  ( 5 ) 

0.120 

0 .086 

t  (6) 

0.120 

0.085 

t  (7) 

0.120 

0.086 

t  { 8 ) 

0.120 

0.086 

t  (9) 

0.120 

0 .086 

t  (10) 

0.120 

0.085 

Vol 

363.916 

170.396 

t : 

plate  thickness  (in) 

Vol : 

volume  (in 

2.  Case 

IB  with  Uniform  Distributed 

Load 

2 

A  uniform  distributed  load,  0.278  lbs/in  was  applied 
over  the  whole  plate.  The  final  optimization  results  show 
that  there  are  6  critical  deflection  constraints  of  the 
nodes  3,  4,  5,  9,  10  and  15,  and  that  the  plate  thickness 
has  a  simila-  trend  to  case  IA  where  diagonal  elements  are 
thicker  than  those  of  the  off-diagonal  elements.  The  final 
volume  is  170.396  in'*.  The  results  are  summarized  in  Table  I 
with  those  of  case  IA. 

B.  DESIGN  CASE  II — THE  STIFFENED  PLATE 

A  quarter  of  the  stiffened  plate  simply  supported,  as 
shown  in  Fig.  4.3,  was  modelled  4x4  mesh  of  the  plate  ele¬ 
ments  and  12  stiffener  elements. 
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1 .  Case  IIA--Rectangular  Type  Frame  Stiffener 

There  are  17  design  variables:  a  plate  thickness, 
and  S  sets  of  stiffener  heights  and  thicknesses,  as  tabu¬ 
lated  in  Table  II.  The  volume  of  the  stiffened  plate  was 
minimized  subject  to  a  total  of  222  constraints:  stress, 
stiffener  height  to  thickness  ratio,  and  nodal  deflection 
constraints. 

a.  Case  IIAl  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  was  applied  at 
the  center  of  the  whole  plate.  There'  are  14  critical  con¬ 
straints:  the  height  to  thickness  ratios,  and  6  deflections 

of  the  nodes  4,  5,  10,  29,  30  and  31.  The  final  volume  is 
102.017  in3  with  a  uniform  plate  thickness  of  0.051  in.  The 
results  are  summarized  in  Taole  III  for  comparing  with  the 
other  cases. 

b.  Case  I IA2  with  Uniform  Distributed  Load 

2 

A  ’uniform  distributed  loud,  0.278  lb/ in  was 
applied  over  the  whole  plate.  There  are  12  critical  con¬ 
straints:  the  height  to  thickness  ratios,  and  4  deflections 

of  the  plate  nodes  7,  8,  14  and  19.  The  final  volume  is 
65.984  in^  with  a  uniform  plate  thickness  of  0.059  in.  The 
results  of  this  case  are  summarized  in  Table  III  with  those 
of  the  case  IIAl.  The  volumes  of  the  above  two  cases  are 
significantly  reduced  in  comparison  to  those  of  the  unstif¬ 
fened  plate. 
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TABLE  II 


Design  Case 

IIA — Design 

Variables  and 

Nodal 

Connectivity 

element 

design 

node 

number 

number 

variables 

Nl 

N2 

N3 

N4 

pi 

t 

1 

6 

7 

2 

?2 

t 

2 

7 

8 

3 

P3 

t 

3 

8 

9 

4 

p4 

t 

4 

9 

10 

5 

p  5 

t 

6 

11 

12 

7 

p6 

♦- 

U 

7 

12 

13 

8 

p7 

t 

8 

13 

14 

9 

p8 

t 

9 

14 

15 

10 

p9 

t 

11 

16 

17 

12 

plO 

t 

12 

17 

18 

13 

pH 

t 

13 

18 

19 

14 

P12 

t 

14 

19 

20 

15 

pl3 

t 

16 

21 

22 

17 

pl4 

t 

17 

22 

23 

18 

pl5 

t 

18 

23 

24 

19 

pl6 

t 

19 

24 

25 

20 

si 

t^l)  ,h(l) 

30 

35 

s2 

V2)  'h(2} 

35 

36 

s3 

tw(3)  ,h(3) 

36 

37 

s4 

tw(4)  ,  h  ( 4  ) 

37 

38 

s5 

tw(5)  , h ( 5 ) 

29 

30 

s6 

^(6)  ,  h  (6  ) 

28 

29 

s7 

tw(7)  ,  h  ( 7 ) 

27 

28 

s8 

t^(8)  ,  h  ( 8 ) 

26 

27 

s9 

^(5)  ,  h  ( 5 ) 

30 

31 

slO 

t^(6) , h{6) 

31 

32 

sll 

^(7)  , h (7) 

32 

33 

s!2 

^(8)  , h ( 8 ) 

33 

34 
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TABLE  III 


Summary  of  Design  Case  IIA 


variable 

concentrated 

distr ibut 

load 

load 

t 

0.051 

0.059 

tw(1) 

0.438 

0.209 

V2) 

0.318 

0.168 

tw(3) 

0.553 

0.267 

V4) 

0.673 

0.344 

tw(5) 

0.097 

0.138 

V6) 

0.096 

0  .164 

tw<7, 

0.  095 

0.171 

CO 

> 

0.095 

0.155 

h  (1) 

2.189 

1.044 

h  (2) 

1.568 

0.832 

h  (3) 

2.753 

1.329 

h  (4) 

3.  343 

1.709 

h  (5) 

0.486 

0.687 

h  (6) 

0.482 

0.821 

h  (7) 

0.478 

0.851 

h  (8) 

0.477 

0.776 

Vol 

102.017 

65.984 

t : 

plate  thickness  (in) 

V 

frame  thickness  (in) 

h: 

frame  height  (in) 

Vol: 

volume  (inJ) 

There  are  33  design  variables:  a  plate  thickness,  8 
sets  of  the  web  heights  and  thicknesses,  and  8  sets  of  the 
flange  widths  and  thicknesses,  as  tabulated  in  Table  IV. 

The  volume  of  the  stiffened  plate  was  minimized  subject  to 
a  total  of  230  constraints:  stress,  web  height  to  thickness 
ratio,  flange  width  to  thickness  ratio,  and  nodal  deflection 
constraints. 

a.  Case  IIB1  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  was  applied  at 
the  center  of  the  whole  plate.  There  are  12  critical  con¬ 
straints:  the  web  height  to  thickness  ratios,  the  3rd  and 

8th  flange  width  to  thickness  ratios,  and  deflection  con¬ 
straints  at  the  node  5  of  the  element  p5  and  node  30  of  the 

3 

element  si.  The  final  volume  is  55.091  in  with  a  uniform 
plate  thickness  of  0.023  in.  The  results  of  this  case  are 
summarized  in  Table  V  for  comparison  w.th  the  other  cases. 

b.  Case  IIB2  with  Uniform  Distributed  Load 

2 

A  uniform  distributed  load,  0.278  lbs/in  was 
applied  over  the  whole  plate.  There  are  16  critical  con¬ 
straints:  the  web  height  to  thickness  ratios  except  the 

8th  one,  the  3rd  through  8th  flange  width  to  thickness  ratios 
except  the  6th  one,  and  4  deflections  of  the  plate  nodes  7, 

3,  14  and  19.  The  final  volume  is  64.414  in^  with  a  uniform 
plate  thickness  of  0.060  in.  This  volume  is  greater  than  the 
volume  of  the  above  concentrated  load  case,  being  different 
from  the  other  cases  which  the  latter  is  greater  than  the 


TABLE  IV 


Design  Case  IIB — Design  Variables  and  Nodal  Connectivity 


element 

design 

node 

number 

number 

variables 

Nl 

N2 

N3 

N4 

Pi 

t 

1 

6 

7 

2 

p2 

t 

2 

7 

8 

3 

P3 

t 

3 

8 

9 

4 

p4 

t 

4 

9 

10 

5 

p3 

t 

6 

11 

12 

7 

p6 

t 

7 

12 

13 

8 

P7 

t 

8 

13 

14 

9 

?3 

t 

9 

14 

15 

10 

P9 

t 

11 

16 

17 

12 

plQ 

t 

12 

17 

18 

13 

pll 

t 

13 

18 

19 

14 

pl  2 

t 

14 

19 

20 

15 

p!3 

t 

16 

21 

22 

17 

pl4 

t 

17 

22 

23 

18 

pi  5 

t 

18 

23 

24 

19 

pi  6 

t 

19 

24 

25 

20 

si 

tw(l) ,h(l)  ,tf (l) , w(l) 

30 

35 

s2 

tw(2) ,h(2) ,tf (1) ,w<l) 

35 

36 

s3 

^(3)  ,h  (  3)  ,tf  (3)  ,  w ( 3) 

36 

37 

s4 

^(4)  ,h(4)  ,tf  (4)  ,w  (4) 

37 

38 

s5 

tw(5) ,h(5) , tf (5) ,w (5) 

29 

30 

so 

tw(6) ,h(6) ,tf (6) , w (6 ) 

28 

29 

si 

tw(7) ,h<7) ,tf (7) ,w(7) 

27 

28 

s8 

tw(8) ,h(8) ,tf (8) ,  w  ( 8) 

26 

27 

s9 

tw(5) ,h(5) ,tf (5) , w ( 5 ) 

30 

31 

slO 

tw(6)  ,h(6)  ,tf  (6)  ,w(6) 

31 

32 

si  1 

t^(7) ,h(7) , tf (7) , w ( 7) 

32 

33 

sl2 

tw(8) ,h(8) ,tf (8) ,w(8) 

33 

34 

44 


TABLE  V 


variable 


tw(D 

tw(2) 

fcW(3) 

V4) 

tT,(5) 

fcW(6) 
tw(7) 
tyj  ( 8) 
h  (1) 
h  (2) 


h 

h 

h 

h 

h 

h 

t 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 
=  (D 


tf  (2) 

tf  (3) 
tf  (4) 
tf  (5) 
tf  (6) 
tf  (7) 

tf  (8) 


Summary  of  Design  Case 

concentrated 

load 

0.023 

0.133 

0.077 

0.137 

0.162 

0.381 

0.355 

0.316 

0.274 

0.659 

0.386 

0.680 

0.805 

1.902 

1.767 

1.570 

1.368 

0.021 

0.028 

0.011 

0.026 

0.876 

0.799 

0.671 

0.460 


IIB 

distributed 

load 

0.060 

0.149 

0.126 

0.197 

0.242 

0.103 

0.114 

0.115 

0.122 

0.736 

0.631 

0.976 

1.201 

0.509 

0.566 

0.569 

0.602 

0.262 

0.183 

0.308 

0.460 

0.209 

0.231 

0.257 

0.201 
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TABLE  V  (Cont'd) 


variable 

concentrated 

load 

distributed 

load 

w  (1) 

0.083 

0.278 

w  (2) 

0.054 

0.187 

w  (3) 

0.057 

0.307 

w  (4) 

0.074 

0.460 

w  (5) 

0.891 

0.210 

w  (6) 

0.814 

0.235 

w  (7) 

0.685 

0.258 

w  (8) 

0.464 

0.201 

Vol 

55.091 

64.414 

t:  place  thickness  (in) 

tTT:  web  thickness  (in) 

w 

h:  web  height  (in) 

tf:  flange  thickness  (in) 

w:  flange  width  (in) 

Vol :  volume  (inJ) 
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former.  This  suggests  that  the  stiffener  type  and  arrange¬ 
ment  used  here  is  particularly  efficient  for  the  concentrated 
load  case.  The  results  of  this  case  are  summarized  in  Table 
V  with  those  of  the  case  IIB1.  The  volumes  of  the  above 
cases  are  significantly  reduced  with  comparison  to  those  of 
unstiffened  plates.  The  volume  of  T  type  stiffener  case 
is  slightly  smaller  than  that  of  rectangular  case  in  the 
distributed  load.  On  the  other  hand,  the  volume  of  T  type 
stiffener  case  is  much  smaller  than  the  rectangular  case 
for  the  concentrated  load. 

C.  DESIGN  CASE  III — THE  HATCH  COVER 

A  quarter  of  a  hatch  cover  was  modelled  4  4  mesh  of  the 

plate  elements  and  24  stiffener  elements,  as  shown  in  Fig. 
4.4. 


1 .  Case  II IA — Rectangular  Type  Frame  Stiffener 

There  are  9  design  variables:  a  plate  thickness, 
and  4  sets  of  stiffener  heights  and  thicknesses,  as  tabulated 
in  Table  VI.  The  volume  of  the  hatch  cover  was  minimized 
subject  to  a  total  of  274  constraints:  stress,  height  to 
thickness  ratio,  and  nodal  deflection  constraints, 
a.  Case  IIIA1  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  was  applied  at 
the  center  of  the  whole  plate.  There  are  4  critical  con¬ 
straints:  the  2nd  and  4th  height  to  thickness  ratios,  and 

the  upper  and  lower  stresses  of  the  node  5  of  the  element 
p4 .  The  dimensions  of  the  2nd  stiffener  set  are  negligibly 
small  with  comparison  to  those  of  the  other  stiffener  sets. 
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TABLE  VI 

Design  Case  IIIA — Design  Variables  and  Nodal  Connectivity 


element 

design 

node 

number 

number 

variables 

Nl 

N2 

N3 

N4 

Pi 

t 

1 

6 

7 

2 

P  2 

t 

2 

7 

8 

3 

p3 

t 

3 

8 

9 

4 

p4 

t 

4 

9 

10 

5 

p5 

t 

6 

11 

12 

7 

p6 

t 

7 

12 

13 

3 

P7 

t 

8 

13 

14 

9 

p8 

t 

9 

14 

15 

10 

p9 

t 

11 

16 

17 

12 

plO 

t 

12 

17 

18 

13 

oil 

t 

13 

18 

19 

14 

P12 

t 

14 

19 

20 

15 

pi  3 

t 

16 

21 

22 

17 

pi  4 

t 

17 

22 

23 

18 

pi  5 

t 

18 

23 

24 

19 

pi  6 

t 

19 

24 

25 

20 

si 

l/2tw(3)  ,  h (  3) 

26 

27 

s2 

l/2tw(3) / h ( 3) 

27 

28 

s3 

tw(l) ,h(l) 

28 

29 

s4 

tw(l) ,h(l) 

29 

30 

s5 

tw(4) , h ( 4 ) 

26 

31 

s6 

tw(2) ,h (2) 

28 

32 

s7 

tw(l) , h ( 1) 

30 

33 

s8 

tw(4)  , h  (4 ) 

31 

34 

s9 

tw(2) ,h(2) 

32 

36 

slO 

tw(l) f h (1) 

33 

38 

si  1 

tw(3)  , h  ( 3) 

34 

35 

sl2 

tw(3)  , h ( 3) 

35 

36 
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TABLE  VI  (Cont'd) 


element 

design 

number 

variables 

Nl 

si  3 

tw(2) ,h (2) 

36 

s!4 

t^(2)  ,h(2) 

37 

si  5 

tw(4)  , h  ( 4 ) 

34 

sl6 

tw(3)  ,h{3) 

36 

si  7 

l/2tw(3)  ,  h  (  3) 

38 

sl8 

tw(4)  /  h  (4 ) 

39 

sl9 

tw(3)  ,h(3) 

40 

s20 

l/2tw(3) , h ( 3) 

41 

s21 

tw(4)  /  h  (4  ) 

42 

s22 

tw(4) , h (4 ) 

43 

s2  3 

tw(4)  ,h(4) 

44 

s24 

tw(4) , h (4 ) 

45 

node  number 
N2  N3  N4 

37 

38 

39 

40 

41 

42 
44 
46 

43 

44 

45 

46 
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This  suggests  that  the  1st  stiffener  set  is  stiff  enough  tc 

compensate  the  elements  attached  to  the  2nd  stiffener  set. 

2 

The  final  volume  is  323.639  in  with  a  uniform  plate  thick¬ 
ness  of  0.359  in.  The  results  of  this  case  are  summarized 
in  Table  VII. 


TABLE  VII 

Summary  of  Design  Case  IIIA 


variable 

concentrated 

load 

distribut 

load 

t 

0.359 

0.059 

twU) 

0 . 811 

0.663 

V2’ 

0.001 

0.002 

V3> 

0.084 

0.518 

0.015 

0.009 

h  (1) 

1.306 

1.750 

h  (2) 

0.006 

0.011 

h  (3) 

0.162 

1.484 

h  (4) 

0.073 

0.019 

Vol 

323.639 

87.314 

t : 

plate 

thickness  (in) 

fcW: 

frame 

thickness  (in) 

h: 

frame 

height  (in) 

Vol: 

volume  (in0) 

b.  Case  IIIA2  with  Uniform  Distributed  Load 

2 

A  uniform  distributed  load,  0.278  lbs/ in  was 
applied  over  the  whole  plate.  There  are  6  critical  constraints 
the  2nd  height  to  thickness  ratio,  and  5  deflections  of  the 
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plate  nodes  8,  9  and  14  and  of  stiffener  nodes  32  and  37. 

With  the  same  suggestion  as  the  concentrated  load  case  the 
dimensions  of  the  2nd  stiffener  set  are  negligibly  small. 

3 

The  final  volume  is  87.314  in  with  a  uniform  plate  thickness 
of  0.059  in.  The  design  results  of  the  above  both  cases 
are  summarized  in  Table  VII  with  those  of  case  IIIA1.  The 
volumes  of  the  above  two  cases  are  significantly  reduced 
with  comparison  to  those  of  unstiffened  plates. 

2 .  Case  IIIB — T  Type  Frame  Stiffener 

There  are  17  design  variables:  a  plate  thickness,  4 
sets  of  the  web  heights  and  thicknesses,  and  4  sets  of  the 
flange  widths  and  thicknesses,  as  tabulated  in  Table  VIII, 

The  volume  of  the  hatch  cover  was  minimized  subject  to  a  total 
of  278  constraints:  stress,  web  height  to  thickness  ratio, 
flange  width  to  thickness  ratio,  and  nodal  deflection 
constraints. 

a.  Case  IIIB1  with  A  Concentrated  Load 

A  concentrated  load  of  1,000  lbs  was  applied  at 
the  center  of  the  whole  plate.  There  are  3  critical  con¬ 
straints:  the  3rd  flange  width  to  thickness  ratios,  and  2 

deflections  of  the  plate  node  5  and  stiffener  node  30. 
Negligibly  small  are  the  whole  dimensions  of  the  2nd  stiffener 
set  and  the  flange  dimensions  of  the  3rd  stiffener  set  with 
comparison  to  the  other  dimensions.  The  final  volume  is 
95.883  in^  with  a  uniform  plate  thickness,  0.083  in. 
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TABLE  VIII 

Design  Case  IIIB — Design  Variables  and  Nodal  Connectivity 


element 

design 

node 

number 

number 

variables 

N1 

N2 

N3 

N4 

Pi 

t 

1 

6 

7 

2 

P2 

t 

2 

7 

8 

3 

p3 

t 

3 

8 

9 

4 

p4 

t 

4 

9 

10 

5 

P5 

t 

6 

11 

12 

•7 

/ 

p6 

t 

7 

12 

13 

8 

t 

8 

13 

14 

9 

p8 

t 

9 

14 

15 

10 

p9 

t 

11 

16 

17 

12 

plO 

t 

12 

17 

18 

13 

pH 

t 

13 

18 

19 

14 

pi  2 

t 

14 

19 

20 

15 

P1 3 

t 

16 

21 

22 

17 

pi  4 

t 

17 

22 

23 

18 

pi  5 

t 

18 

23 

24 

19 

pi  6 

t 

19 

24 

25 

20 

si 

l/2tw(3)  ,h(3)  ,tf  (3)  , l/2w  ( 3} 

26 

27 

s2 

l/2tw(3)  , h  (  3)  ,tf (3)  , l/2w ( 3) 

27 

28 

s3 

tw ( 1 ) ,h (1) , tf (1) ,w(l) 

28 

29 

s4 

tw(l) , h (1 ) ,tf (1) ,w(l) 

29 

30 

s5 

tw ( 4 ) ,h  ( 4 ) ,tf (4) ,w(4) 

26 

31 

s6 

t^(2)  ,h(2)  ,tf  (2)  ,w(2) 

28 

32 

s7 

tw(l) ,h (1) ,tf (1) ,w(l) 

30 

33 

s8 

tw(4) , h (4 ) ,tf (4) ,w(4 ) 

31 

34 

s9 

tw(2) ,h(2) ,tf (2) ,  w  <  2 ) 

32 

36 

slO 

tw(l) ,h(l) ,tf (1) ,w(l) 

33 

38 

sll 

tw(3)  ,h(3>  ,tf (3)  , w(  3) 

34 

35 

si  2 

tw(3)  ,h(3)  ,tf  (3)  ,w(3) 

35 

36 
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TABLE  VIII  (Cont'd) 


element 

design 

node  number 

numoer 

variables 

Nl 

N2  N3 

sl3 

tw(2) , h { 2) ,tf (2) , w (2) 

36 

37 

sl4 

tw(2) ,h (2) ,tf (2) , w ( 2 ) 

37 

38 

sl5 

tw(4) ,h (4) , tf (4) , w ( 4 ) 

34 

39 

si  6 

tw(3)  ,h  (3)  ,tf  (3)  ,  w  (  3) 

36 

40 

sl7 

l/2tw(3) ,  h  ( 3) ,tf (3) , l/2w ( 3) 

38 

41 

sl8 

tw(4) ,h (4) ,tf (4) , w ( 4 ) 

39 

42 

sl9 

tw(3)  ,  h  (  3)  ,tf  (3)  ,  w  (  3) 

40 

44 

s20 

l/2t^(  3)  ,  h  (  3)  ,tf  (3)  ,  l/2w  (  3) 

41 

46 

s21 

tw(4) ,  h  ( 4  ) ,tf (4) ,w (4) 

42 

43 

s22 

tw(4) ,h(4) ,tf (4) , w (4 ) 

43 

44 

s2  3 

tw(4)  ,h(4)  ,tf  (4)  ,  w  ( 4  ) 

44 

45 

s24 

tw(4) ,h (4)  ,tf  (4) , w ( 4 ) 

45 

46 
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b.  Case  IIIB2  with  Uniform  Distributed  Load 

2 

A  uniform  distributed  load,  0.278  lbs/in  was 
applied  over  the  whole  plate.  There  are  9  critical  con¬ 
straints:  the  3rd  web  height  to  thickness  ratio,  the  3rd 

flange  width  to  thickness  ratio,  4  deflections  of  the  plate 
nodes  8,  9,  13  and  14,  and  3  deflections  of  the  stiffener 
nodes  32,  36  and  37.  Negligibly  small  are  the  whole  dimen¬ 
sions  of  the  2nd  stiffener  set,  and  the  flange  dimensions  of 
the  3rd  and  4th  stiffener  sets  with  comparison  to  the  other 
dimensions.  This  suggests  that  the  1st  stiffener  set  is 
stiff  enough  to  compensate  the  elements  attached  to  the  2nd 
stiffener  set.  Additionally  it  can  be  suggested  that  the 
3rd  and  4th  stiffener  sets  are  stiff  enough  to  support  the 
given  load  with  only  the  rectangular  type  ones.  The  final 
volume  is  76.436  in'3  with  a  uniform  plate  thickness,  0.062 
in.  The  results  of  the  above  both  cases  are  summarized  in 
Table  IX.  The  volumes  of  the  above  two  cases  are  signifi¬ 
cantly  reduced  with  comparison  to  those  of  the  unstiffened 
plates  and  the  volume  of  T  type  stiffener  case  is  much 
smaller  than  that  of  the  rectangular  type  case. 


55 


■■ 


TABLE  IX 


Summary  of  Design  Case  IIIB 


variable 

concentrated 

load 

t 

0.083 

Vi} 

0.284 

v2) 

0.003 

V3> 

0.323 

V4> 

0.020 

h  (1) 

1.201 

h  (2) 

0.010 

h  (3) 

1.433 

h  (4) 

0.060 

tf  (1) 

1.001 

tj  (2) 

0.003 

t£  (3) 

0.000 

t  -  ( 4) 

0.010 

w‘(l) 

1.110 

w  (2) 

0.010 

w  (3) 

0.002 

w  (4) 

0.0  30 

Vol 

95.883 

t : 

plate  thickness  (in) 

^ : 

web  thickness  (in) 

h : 

web  height  (in) 

tV 

flange  thickness  (in 

w: 

flange  width  (in) 

Vol: 

volume  (in2) 

distributed 

load 

0.062 

0.284 

0.003 

0.318 

0.020 

1.200 

0.010 

1.430 

0.060 

1.000 

0.003 

0.001 

0.010 

1.110 

0.010 

0.003 

0.001 

76.436 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

A  finite  element  analysis  program  was  developed  and 
coupled  to  an  optimizer  for  stiffened  plate  design.  Stif¬ 
fened  plates  were  designed  for  minimum  volume  with  a  concen¬ 
trated  load. 

Five  important  conclusions  are  made  from  considering 
several  design  examples. 

First,  a  stiffened  plate  is  much  more  efficient  than 
the  unstiffened  plate  for  the  optimum  volume  design. 

Second,  it  is  required  to  attach  the  larger  stiffeners 
to  the  diagonal  elements  of  the  plate  as  compared  to  the 
off-diagonal  elements. 

Third,  the  T  type  stiffeners  provide  a  lower  volume 
design  than  the  rectangular  type  stiffeners. 

Fourth,  the  stiffened  plate  with  the  diagonally  attached 
stiffeners  are  much  more  efficient  than  those  with  the 
rectangularly  attached  stiffeners. 

Fifth,  the  stiffened  plate  attached  diagonally  with  the 
T  type  stiffeners  is  particularly  efficient  for  the  concen¬ 
trated  load  case. 

These  conclusions  are  based  on  the  design  of  simply 
supported  plates  and  the  conclusions  could  be  quite  differ¬ 
ent  if  other  boundary  conditions  are  considered. 
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B.  RECOMMENDATIONS 


The  following  recommendations  may  be  useful  for  future 
work . 

1.  The  analysis  program  should  be  modified  to  a 
generalized  program  with  a  variety  of  elements. 

2.  The  analysis  program  should  be  extended  to 
consider  the  dynamic  problem. 

3.  The  analysis  program  should  be  extended  to 
consider  the  multiple  loading  conditions  . 

4.  Routines  should  be  added  to  calculate  gradients 
analytically . 

5.  Buckling  constraints  of  the  plate  and  stiffener 
elements  should  be  added  in  the  optimization 
process . 

6.  Frequency  constraints  should  be  added  in  the 
optimization  process. 

7.  The  analysis  program  should  be  modified  to  reduce 
the  computing  time  during  the  optimization  process. 
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